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We obtain a new exact equilibrium solution to the TV-body problem in a one-dimensional rela- 
tivistic self-gravitating system. It corresponds to an expanding/contracting spacetime of a circle 
with TV bodies at equal proper separations from one another around the circle. Our methods 
are straightforwardly generalizable to other dilatonic theories of gravity, and provide a new class of 
solutions to further the study of (relativistic) one-dimensional self-gravitating systems. 
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Circular gravity is a general relativistic gravitational system in (1 + 1) dimensions in which the topology of space 
is equivalent to that of a circle S 1 . Its counterpart, lineal gravity (the same system with the S* 1 topology replaced 
by M 1 ), has been the subject of much investigation over the past twenty years, motivated by the desire to understand 
black holes [jjj , spacetime structure [^],[| , canonical quantum gravity [Q , string physics || , information loss || and the 
TV-body problem [0J|]. This latter problem involves consideration of a collection of TV particles interacting through 
their own mutual gravitational attraction, along with other specified forces. Apart from functioning as prototypes 
for higher dimensional gravity, in lineal gravity such self-gravitating systems also approximate the behaviour of some 
. | physical systems in 3 spatial dimensions, such as stellar dynamics orthogonal to a galactic plane, the collisions of flat 
i parallel domain walls, and the dynamics of cosmic strings. The ergodic and equipartition properties of non- relativistic 
one-dimensional self-gravita ting systems are a current subject of research j|, and recently such systems were shown 
fSJ , to exhibit fractal behaviour |H]] . In the past few years these systems have been extended to relativistic lineal gravity. 
£> ■ The 2-body problem has been solved exactly JO, and has been extended to include both cosmological expansion 
0^ | |£2|,HI and electromagnetic interactions JT3. The TV-body extension of these systems, while considerably more 
complicated, is nevertheless amenable to analysis [|| and includes the non-relativistic self-gravitating systems [p||l5|] 
as a limiting case. 

Comparitivcly little is known about the circular counterparts of these systems; they introduce qualitatively new 
' features that are absent in the non-compact lineal case. For example, in the absence of other fields there are no 
static non-relativistic equilibrium solutions for TV particles on a circle. This is because the non-relativistic potential 
grows linearly with increasing distance from the source, and a circular topology admits no solutions to the matching 
qh| conditions for any TV > 1. 

We present here an exact equilibrium solution to the TV-body problem in circular gravity. Our solution corresponds 
^jrj to TV equal mass particles at equidistant proper separations from one another in a spacetime that expands/contracts in 
response to these sources. The configuration we find therefore corresponds to a dynamic, albeit unstable, equilibrium. 
To our knowledge these are the first TV-body dynamic equilibrium solutions in any relativistic theory of gravity. While 
our solutions depend upon the choice of dilaton potential, our approach is straightforwardly generalizable to any 
dilatonic theory of gravity, and provides a new class of solutions to further the study of (relativistic) one-dimensional 
5^ ! self-gravitating systems. 

The action integral for a generic dilaton theory of gravity |^],[l6| minimally coupled to TV point masses is 
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where VP is the dilaton field, g^ is the metric with determinant g and covariant derivative V M , R is the Ricci scalar, 
and r a is the proper time of a-th particle, with k = 87rG/c 4 . For circular gravity the range of the spatial coordinate x 
is — L < x < L, and the circular topology implies that all fields must be smooth (or at least C 1 ) functions of x with 
period 2L, which implies 

f(L) = /(-£) and f'(L) = /'(-£) = . (2) 

for all functions. Henceforth we set V (W) = A/2 throughout, where A is constant. 

The action (||) describes a generally covariant self-gravitating system (without collisional terms, so that the bodies 
pass through each other); its field equations are 
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where the stress-energy due to the point masses is 

T /^ = X] ma j dT a^=^9^a9up'^- l ^-5 2 (x - Z a (T a )) (5) 

and is conserved. The preceding system generalizes Jackiw-Teitelboim lineal gravity |jl8| , which equates scalar curva- 
ture to a (cosmological) constant A, reducing to it if the stress-energy vanishes. The specific choice of dilaton potential 
ensures that the evolution of the dilaton does not modify the aformentioned reciprocal gravity/matter dynamics. The 
system (|J) is a closed system of N + 1 that describes the evolution of the single metric degree of freedom and the N 
degrees of freedom of the point masses. The evolution of the point-masses governs the evolution of the dilaton via 
([|) , whose divergencelessness is consistent with the conservation of T^ v , yielding only one independent equation to 
determine the single degree of freedom of the dilaton. 

Making use of the previously formulated general framework for canonical circular dilaton gravity Jl7| , we write the 

metric as ds 2 = — (Nodt) 2 + 7 (dx + ^dtj ■ The extrinsic curvature K is then 

K = (2iV o7 )- 1 (2a 1 7V 1 - 7- 1 JVi0i7 - d al ) = ^k(tt - II/7) (6) 

where the conjugate momenta to 7 and 'J are respectively denoted 7r and II. We can take the extrinsic curvature 
if to be a time coordinate r (t) of the circle, thereby allowing the elimination of tt from all field equations and the 
identification of the Hamiltonian from the action (Q) after canonical reduction jlTj : 



I = I d-.r { > j,„ :„<M.r - z a ) + II-| + In 7) - H I , (7) 



where H = J dxTL =— J dx^/j, which is the circumference functional of the circle when f is constant. The dynamics 
is that of a time-dependent system, with the time-dependence corresponding to the time-varying circumference of 
the circle of constant mean extrinsic curvature. As with general relativity on spatially compact manifolds in (2 + 1) 



dimensions |19|, the spatial metric can be chosen so that 7 = j(t) , i.e. 7' = , since any metric on a circle is 
globally conformal to a flat metric. Choosing a time parametrization so that is constant, the Hamiltonian will 
be time- independent. 



We find after some manipulation that the canonical form of the field equations (|3|,[|) is 

211' - IW' = 0; (8) 
¥" - i(*') 2 - («n) 2 + 7 (r 2 - A/2) + k Vl m 5{x - z a (x")) = (9) 

a 

N{ -t^-NoI^ -A/2)i + K /2Y^/m5(x-z a (x ))\ =0 (10) 



N[ — 7/2 + 7TA/0 = (11) 

tl + d 1 (--N 1 U+-^—N ^ + ^—N^ = (12) 
7 2k^/7 k^/7 

i> + 2N ( k^ + t ) -JVi(-*') = (13) 

where equilibrium solutions are characterized by z a = = p a , corresponding to a situation in which all particles are 
motionless at various points around the circle. A calcuation of the geodesic equations in (||) show that this implies 

^ =0 ATi(z Q )=0 . (14) 
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Remarkably, this set of equations admits an exact solution for N bodies of equal mass m at equal time- varying proper 
separations. We set c 2 = 7 (t 2 - A/2) = c? + > 0, and r = \z a — z +i| = ^ > where a = 1,2, . . . ,N. The system is 
cyclically symmetric and so we can choose the origin to be halfway between any two particles in the N =even case or on 
a particle in the N =odd case. In seeking solutions to eqs. ds|-|lS|) we must ensure that the solutions for \& and No are 
appropriately matched at the locations z a of each particle and the identification point |x| = L] they must be continuous 
there but not differentiable at the locations of the particles; for example lim e ^o \^'{z a + e) — ^'(z a — e)] = n^/jm. 

With this in mind, eq. (||) is straightforwardly found to yield II = Tlo(t) exp where Ilo(t) is an arbitrary 
function. Upon substitution of this result into eq. @ , and employing the conditions (|lj) , the next three equations 
are solved by straightforward integration and matching across the particles, furnishing a solution for "f, No, and Ni 
in terms of n (t) and the metric function 7 (t). Imposing ( |l4| ) and the periodicity condition N-^(L) = Nx(—L) then 
yields a differential equation for 7 whose solution is 

LsJt? - A/2 \ <gff + P(T 2 -A/2) ' 1 



where M = mN is the total mass of the system. The function IIo (t) is then found from integration of ([13|). Satisfaction 
of equation ( |l2| ) provides a non-trivial consistency check on the full result, which is 

* = - 21 " (» ~ iss^p*^ - -1 - L ») - 2 ^7^t ) (16) 

N " = ^ (sh§S - -\ - *>» - ") (") 

N X = J^x - 7 3 7 T Sm , )~> z ^ i^Hc+(\x z a \ L)) + sinh(c + L)] (18) 
c_|_ Cj^Ij siH-ii^c_|_xvj _^ 

for the dilaton, lapse and shift functions, where the step function e (x) = — , and e (0) = 0. Here 
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no = ^( T ^ 16A ,f ( y_ 1) + (^-A/2)j (20) 
with £ and £ integration constants. The conjugate dilaton momentum is easily seen to be 



*n= _, (21, 

s3=5}E:„«».h(c + (ii-^i-L))-/3 

We can more succinctly write the above using the relation cosh(c+(|x — z a \ — L)) = cosh(c + /(a;)) smh ^+-^) 

sinh(^t_) 

where f(x) is the saw-tooth function that peaks with a value of L/N (i.e. f(z a ) — L/N) at the particle locations and 
has zero magnitude half-way between the particles (i.e. j( z °~^°+ 1 ) = ). This relation holds because each function 
is periodic, and so the origin can be shifted to be between any pair of adjacent particles. A simple shift of the origin 
in the cosh function and a subsequent manipulation of the sum yields the equivalence. 

Another class of solutions exists for c 2 = c 2 _ < 0. These are obtained upon replacing c + — > ic_ in the above results. 
The solution for the metric function 7 now becomes 



r N 

V7 = 



L^/K/2-t 2 



^ + (l-f)(A/2-rVf 



arctan — =-77= \/A/2 — r 2 + kn 

1 ^_£2 (A/2 _ r2) V ! 



(22) 



where k is an integer and the integration constant |£| < y 1 + g^g^s as a consequence of the periodicity of Ni(x,t). 
There are no solutions for c 2 = 0. 
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The solution (fL5|- |l8| , [2l| ) (and its c < counterpart) is a new dynamical equilibrium solution to the iV-body 
problem in relativistic gravity. It corresponds to an expanding/contracting spacetime of a circle with N bodies at 
equal time-varying proper separations from one another around the circle. For N = 1 the results reduce to the single 
particle case obtained previously |fl7| ; indeed, the TV-body solution for the spatial metric is equivalent to that of 
the single-particle solution with a rescaling of the parameters L — ► L/N and M — > M/N. Note that this rescaling 
equivalence does not hold for the remaining metric and dilaton functions. 

Analysis of the behaviour of the spatial metric and Hamiltonian thus proceeds along the same lines as ref. jl7| . For 
fixed particle mass m the spatial metric and Hamiltonian behave analogously to the single-particle case times a factor 
of N. For A < the proper circumference of the circle expands from zero to some maximal size and then recontracts, 
whereas for A = it undergoes a perturally decelerating expansion due to the presence of the point masses. For 
A > the cosmological expansion opposes decelerating effects due to the point masses and a variety of behaviours 
can ensue, including expansions from zero size that asymptote to some finite value, evolutions to maximal size that 
begin/end at finite circumference and perpetually oscillating universes. 

For fixed total mass M the evolution of the spatial metric becomes independent of M as the number of bodies N 
becomes large, and the full solution more closely approaches its vacuum state on a circle with L — ► L/N. 

Generalization of our approach to other dilatonic theories of gravity is straightforward. An arbitrary (1 + 1)- 
dimensional dilatonic theory of gravity minimally coupled to N bodies can, after conformal transformations and a 
rescaling of the dilaton, be written in the form ([!]) |i],|l6| with the replacement m — > m(ip). The field equations remain 



will modify 



the same, except that ( p"^ ) will contain terms proportional to Noj^ and TVo^S and 

(|l4|). While solving the field equations between the particles is straightforward (20), implementation of the matching 
conditions is not, and so the existence of an equilibrium solution in the general case remains an interesting open 
question. 

The solution we have obtained depends quite sensitively on the matching conditions, and these in turn depend upon 
the equality of the point masses and their proper separations. Perturbations from this equilbrium case would be 
interesting to study, as they would form model inhomogeneous self-gravitating cosmological systems. For sufficiently 
large positive A, we expect that cosmological expansion will overcome the tendency to gravitationally clump, with the 
bodies maintaining some form of thermal equilibrium diluted by the overall expansion. For negative A, there should 
be an instability to formation of a 2d black hole 
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